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ABSTRACT 

Context. The CoRoT satellite has recently discovered the transits of a telluric planet across the disc of a late-type magnetically active 
star dubbed CoRoT-7, while a second planet has been detected after filtering out the radial velocity (hereafter RV) variations due to 
stellar activity. 

Aims. We investigate the magnetic activity of CoRoT-7 and use the results for a better understanding of the impact of magnetic 
activity on stellar RV variations. 

Methods. We derive the longitudinal distribution of active regions on CoRoT-7 from a maximum entropy spot model of the CoRoT 
light curve. Assuming that each active region consists of dark spots and bright faculae in a fixed proportion, we synthesize the expected 
RV variations. 

Results. Active regions are mainly located at three active longitudes which appear to migrate at different rates, probably as a 
consequence of surface differential rotation, for which a lower limit of AQ/C1 = 0.058 ± 0.017 is found. The synthesized activity- 
induced RV variations reproduce the amplitude of the observed RV curve and are used to study the impact of stellar activity on 
planetary detection. 

Conclusions. In spite of the non-simultaneous CoRoT and HARPS observations, our study confirms the validity of the method 
previously adopted to filter out RV variations induced by stellar activity. We find a false-alarm probability < 10~ 4 that the RV 
oscillations attributed to CoRoT-7b and CoRoT-7c are spurious effects of noise and activity. Additionally, our model suggests that 
other periodicities found in the observed RV curve of CoRoT-7 could be explained by active regions whose visibility is modulated by 
a differential stellar rotation with periods ranging from 23.6 to 27.6 days. 

Key words, stars: magnetic fields - stars: late-type - stars: activity - stars: rotation - planetary systems - stars: individual (CoRoT-7) 

1 . Introduction CoRoT-7 is an active star showing a rotational modulation of 

its optical flux wit h an amplitu de up to ~ 0.03 mag and a period 
CoRoT (Convection, Rotation and Transits) is a photometric of ~ 23.5 days (lOueloz et al.ll2009l) . The discovery of CoRoT- 
space experiment devoted to asteroseismology an d the search 7 C and the measurement of the mass of CoRoT-7b have become 
for e xtrasolar planets by the method of transits dBaglin et all possible with the HARPS spectrograph after filtering out the ap- 
\ I2006I) . It has recently discovered two telluric planets around a parent variations of the stellar RV induced by magnetic activity 



X 

j^ . late-type star called CoRoT-7, the innermost of which transits which have an amplitude of » 40 m s , i.e., 5-6 times greater 
across the disc of the starandhas an orbital period of onlyO.8536 than the wobbles produced by the gravitational pull of the plan- 
days, while t he other has a period of 3.69 days an d does not ets. Therefore, a more detailed model of the stellar activity is 
show transits dLeger et al .l l2009t lOueloz et a ll I2009I) . The mass needed to improve RV measurements and to obtain better plan- 
of the inner planet CoRoT-7b is 4.8 ± 0.8 M e , while that of etary parameters. This would also rule out the possibility that a 
CoRoT-7c is 8.4 ± .9 M®, assuming they are on coplanar orbits quasi-periodic RV signal is due to rotating spots on the stellar 
(lOueloz et alj|2009t) . The relative depth of the transit of CoRoT- surface instead of a planetary co mpanion as, e.g., in the case of 
7b is AF/F = (3.35 + 0. 12) x 10~ 4 in the CoRoT white passband HD 166435 (lOueloz et alj|2001h . 
(see Sect. 13, leading to a planetary radius of 1.68 ± 0.09 R®. 

We apply the same spot modelling approach used for the 

high-preci sion light curves of C oRoT-2 dLanza et alj|2009al) and 
* Based on observations obtained with CoRoT, a space project op- CoRoT-4 (|Lanza et al.| |2009bJ), deriving the longitudinal dis- 
erated by the French Space Agency, CNES, with partecipation of the tribution of photospheric brightness inhomogeneities and then- 
Science Programme of ESA, ESTEC/RSSD, Austria, Belgium, Brazil, evolution during the ~ 130 days of CoRoT observations. The 
Germany, and Spain. different rotation rates of active regions allow us to estimate 
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a lower limit for the amplitude of the differential rotation 
of CoRoT-7. From our spot maps we synthesize the activity- 
induced RV perturbations using a simple model for the distortion 
of the line profiles, as described in Sect. 13.21 A Fourier analysis 
of such a synthetic time series yields typical periods induced by 
stellar activity, allowing us to discuss its impact on the detection 
of planets around CoRoT-7. 

2. Observations 

CoRoT-7 has been observed during the first long run of CoRoT 
toward the galactic anticentre from 24 October 2007 to 3 March 
2008. Since the star is bright (V — 11.67), the time sampling 
is 32 s from the beginning of the observations. CoRo T performs 
aperture photometry with a fixed mask (see Fig. 6 in iLeger et all 
2009). The contaminating flux from background stars falling in- 
side the mask amounts to a maximum of 0.9 percent and pro- 
duces a dilution of the photometric variation of CoRoT-7 lower 
than 1.8 x 10~ 4 mag th at can be safely ignored for our purposes 
(see iLeger et aLll2009l for upper limits on the background flux 
variations). The flux inside the star's mask is split along detector 
column boundaries into broad-band red, green, and blue chan- 
nels. 

The ob s ervatio ns and data processing are described by 
ILeger et all (120091) . to whom we refer the reader for details. 
The reduction pipeline applies corrections for the sky back- 
ground and the pointing jitter of the satellite, which is par- 
ticularly relevant during ingress and egress from the Earth 
shadow. Measurements during the crossing of the South Atlantic 
Anomaly of the Earth's magnetic field, which amounts to about 
7-9 percent of each satellite orbit, are discarded. High- 
frequency spikes due to cosmic ray hits and outliers are removed 
by applying a 7-point running mean. The final duty cycle of 
the 32-s observations is 88.8 percent for the so-called N2 data 
time series that are accessible through the CoRoT Public Data 
Archive at IAS (http://idoc-corot.ias.u-psud.fr/). More informa- 
tion on the instrument, its operation, and performance can be 
found in lAuvergne e t al. (2009). 

To increase the signal-to-noise ratio and reduce systematic 
drifts present in individual channels, we sum up the flux in the 
red, green, and blue channels to obtain a light curve in a spectral 
range extending from 300 to 1 100 nm. 

Th e transits are removed using the ephemeris of ILeger et alj 
(2009) and the out-of- transit data are binned by computing 
average flux values along each orbital period of the satellite 
(6184 s). This has the advantage of removing tiny systematic 
variations associated with the orbital motion of the satellite (cf. 
lAlonso et al.ll2008l : lAuvergne et alJl2009h . In such a way we ob- 
tain a light curve consisting of 1793 average points ranging from 
HJD 2454398.0719 to HJD 2454528.8877, i.e., with a duration 
of 130.8152 days. The average standard deviation of the points 
is 2.678 x 10~ 4 in relative flux units. The maximum observed 
flux in the average point time series at HJD 2454428.4232 is 
adopted as a reference unit level corresponding to the unspotted 
star, since the true value of the unspotted flux is unknown. 

3. Light curve and radial velocity modelling 

3.1. Spot modelling of wide-band light curves 

The reconstruction of the surface brightness distribution from 
the rotational modulation of the stellar flux is an ill-posed prob- 
lem, because the variation of the flux vs. rotational phase con- 
tains only information on the distribution of the brightness inho- 



mogeneities vs. longitude. The integration over the stellar disc 
effectively cancels any latitudinal information, particularly when 
the inclination of the rotation axis along the line of sight is close 
to 90°, as it is ass umed in the present case (see Sect. [4] and 
lLanza et aLll20 09a). Therefore, we need to include a priori in- 
formation in the light curve inversion process to obtain a unique 
and stable map. This is done by computing a maximum entropy 
(hereafter ME) map, which has been proven to successfully re- 
produce active region distribution and area variations in the case 
of the Sun (cf. lLanza et al.H"2007l) . 

In our model, the star is subdivided into 200 surface ele- 
ments, namely 200 squares of side 18°, with each element con- 
taining unperturbed photosphere, dark spots, and facular areas. 
The fraction of an element covered by dark spots is indicated by 
the filling factor /, the fractional area of the faculae is Qf, and 
the fractional area of the unperturbed photosphere is 1 -(Q+ l)f. 
The contribution to the stellar flux coming from the k-th surface 
element at the time tj, where jf = 1, ...,N, is an index numbering 
the N points along the light curve, is given by: 



AF kj = /oGu*/>U-(G + l)/ + c 8 /+ 



(1) 



where Iq is the specific intensity in the continuum of the unper- 
turbed photosphere at the isophotal wavelength of the observa- 
tions, c s and Cf are the spot and facular contrasts, respectively (cf. 
lLanza et al]|2004l) . A k is the area of the k-th surface element, 



wQikj) = 



1 ifn kj >0 
if Hkj < 



(2) 



is its visibility, and 



Hkj = cos if/kj = sin i sin 9k cos[4 + Q(f j - to)] + cos i cos 9k, (3) 

is the cosine of the angle t/^ 7 between the normal to the sur- 
face element and the direction of the observer, with i being the 
inclination of the stellar rotation axis along the line of sight, 
6k the colatitude and Ik the longitude of the fc-th surface el- 
ement; Q denotes the angular velocity of rotation of the star 
(Q. = 2n/P ml , with P rot the stellar rotation period), and to the ini- 
tial time. The specific intensity in the continuum varies accord- 
ing to a quadratic limb-darkening law, as adopted by Lanza et al. 
d2003l) for the case of the Sun, viz. Iq °c a p + b p fi + c p fj. 2 . The 
stellar flux computed at the time tj is then: F{tj) - Y,k &Fkj. To 
warrant a relative precision of the order of 10~ 5 in the compu- 
tation of the flux F, each surface element is further subdivided 
into 1° x l°-elements and their contributions, calculated accord- 
ing to Eq. (Q]), are summed up at each given time to compute 
the contribution of the 18° x 18° surface element to which they 
belong. 

We fit the light curve by varying the value of / over the sur- 
face of the star, while Q is held constant. Even fixing the rotation 
period, the inclin ation, and the spot and facular contrasts (see 
lLanza et aDl2007l for details), the model has 200 free param- 
eters and suffers from non-uniqueness and instability. To find 
a unique an d stable spot map, we apply ME regularization, as 
described in lLanza et al] (l2007t) . by minimizing a functional Z, 
which is a linear combination of the y 1 an d the entropy func- 
tional S; i.e., 



Z= X \f)-riS{f), 



(4) 



where / is the vector of the filling factors of the surface ele- 
ments, rj > a Lagrangian multiplier determining the trade-off 
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between light curve fitting and reg ularization; the expression for 
S is given in lLanza et all ([1998). The entropy functional S is 
constructed in such a way that it attains its maximum value when 
the star is free of spots. Therefore, by increasing the Lagrangian 
multiplier 77, we increase the weight of S in the model and the 
area of the spots is progressively reduced. This gives rise to sys- 
tematically negative residuals between the observations and the 
best-fit model when 77 > 0. The optimal value of 77 depends on 
the information content of the light curve, which in turn depends 
on the ratio of the amplitude of its rotational modulation to the 
average standard deviation of its points. In the case of CoRoT-7, 
the amplitude of the rotational modulation is ~ 0.018, while the 
average standard deviation of the points is ~ 2.7 x 10~ 4 in rela- 
tive flux units, giving a signal-to-noise ratio of ~ 65 - 70. This 
is sufficient to adopt the same regularization criterion applied in 
the case of CoRoT-2 (see Lanza et al. 2009a), iterating the value 
of 77 until the condition |/j reg | =* eo is met, where /i leg is the mean 

of the residuals and eo = <t () / y/N their standard error, i.e., the 
ratio of their standard deviation ctq in the case of the unregular- 
ized best fit (i.e., for 77 = 0) to the square root of the number N of 
points in each fitted subset of the light curve having a duration 
Aff (see below). 



In the case of the Sun, by assuming a fixed distribution of the 
filling factor, it is possible to obtain a good fit of the irradiance 
changes only for a limited time interval Aff, not exceeding 14 
days, which is the lifetime of the largest sunspot groups domi- 
nating the irradiance variation. In the case of other active stars, 
the value of Atf must be determined from the observations them- 
selves, looking for the longest time interval that allows a good 
fit with the applied model (see Sect.|Ui. 

The optimal values of the spot and facular contrasts and of 
the facular-to-spotted area ratio Q in stellar active regions are 
unknown a priori. In our model the facular contrast Cf and the 
parameter Q enter as the product C{ Q, so we can fix Cf and vary 
Q, estimating its best value by minimizing the x 1 of the model, 
as shown in Sect. [4] Since the number of free parameters of the 
ME model is large, for this speci fic application we make use of 
the model of lLanza et alj d2003l) . which fits the light curve by 
assuming only three active regions to model the rotational mod- 
ulation of the flux plus a uniformly distributed background to 
account for the variations of the mean light level. This proce- 
dure is the same that was used to fix the value of in the cases 
of CoRoT-2 and CoRoT-4 (cf. lLanza et alJl2009allbh . 



We shall assume an inclination of the rotation axis of C0R0T- 
7 of 7 =s 80° (see Sect. 0). Since the information on spot lati- 
tude that can be extracted from the rotational modulation of the 
flux for such a high inclination is negligible, the ME regulariza- 
tion virtually puts all the spots at the sub-observer latitude (i.e., 
90° - i ~ 10°) to minimize their area and maximize the entropy. 
Therefore, we are limited to mapping only the distribution of the 
active regions vs. longitude, which can be done with a resolution 
better than ~ 50° (cf. Lanza et al. 2007, 2009a). Our ignorance 
of the true facular contribution may lead to systematic errors in 
the active region longitudes derived by our model because fac- 
ulae produce an increase of the flux when they are close to the 
limb, leading to a systematic shift of the longitudes of the ac- 
tive regions used to reproduce the observed flux modulation, as 
discussed by Lanza et al. (2007) for the case of the Sun and il- 
lustrated byJLanzaeFalJ (12009al, cf. Figs. 4 and 5) for CoRoT-2. 



3.2. Activity-induced radial velocity variations 

Surface magnetic fields in late-type stars produce brightness 
and convection inhomogeneities that shift and distort their spec- 
tral line profiles leading to apparent RV variations (cf., e.g., 
ISaar & Donahuelll997t ISaar et :alJH998L iHuber et al.ll2009h . To 
compute the apparent RV variations induced by stellar active re- 
gions, we adopt a simple model for each line profile considering 
the residual profile R(A) at a wavelength A along the line, i.e.: 
R(A) = [1 — I{A)/I C ], where I{A) is the specific intensity at wave- 
length A and I c the intensity in the continuum adjacent to the line. 
The local residual profile is assumed to have a Gaussian shape 
with thermal and microturbulent width A/Id (cf., e.g., iGravl 
1988, Ch. 1), i.e.: 



R(A) oc exp 



A-A 
AA n 



(5) 



where Aq is the local central wavelength. For simplicity, R(A) is 
assumed to be independent of the disc position and the presence 
of spots or faculae, so the local specific intensity along the line 
1(A) depends only on the variation of the continuum intensity 7 C 
owing to limb-darkening and the effects of dark spots and bright 
faculae, as specified by Eq. ([TJ. 

To include the effects of surface magnetic fields on convec- 
tive motions, we consider the decrease of macroturbulence ve- 
locity and the reduction of convective blueshifts in active re- 
gions. 

Specifically, we assu me a local r adial-tangential macroturbu- 
lence, as introduced by Gray dl988l) . with a distribution function 
for the radial velocity v of the form: 



exp 



®(v,n) 



exp 



\ &Tj" 7 



2 V^RT V 1 "A* 2 



V^Trt/j 



(6) 



where £rt is the macroturbulence and fi is given by Eq. (|3). 
The profile emerging from each surface element is obtained by 
convolving the local profile in Eq. ([5]) with the Doppler shift 
distribution as generated by the macroturbulence function in 
Eq. ©. The parameter ^ is assumed to be reduced in spot- 
ted and facular areas, according to their total filling factor, i.e., 
^rt = £6[1 ~~ (Q + !)/]> where £0 is the unperturbed macroturbu- 
lence. This is a convenient parameterization of the quenching ef- 
fect of surface magnetic fields on convective turbulence, at least 
in the framework of our simplified model, and is supported by 
the observations of the reduction of tu rbulent velocity fields in 
the plage regions of the Sun (cf., e.g., iTitle et afl l992). 

Convective blueshifts arise from the fact that in stellar photo- 
spheres most of the flux comes from the extended updrafts at the 
centre of convective granules. At the centre of the stellar disc, the 
vertical component of the convective velocity produces a max- 
imum blueshift, while the effect vanishes at the limb where the 
projected velocity is zero. The cores of the deepest spectral lines 
form in the upper layers of the photosphere where the vertical 
convective velocity is low, while the cores of shallow lines form 
in deeper layers with a greater vertical velocity. Therefore, the 
cores of shallow lin es are blue sh ifted with respect to the cores 
of the deepest lines. IGravl (120091) showed this effect by plotting 
the endpoints of line bisectors of shallow and deep lines on the 
same velocity scale. He showed that the amplitude of the rela- 
tive blueshifts scales with the spectral type and the luminosity 
class of the star. For the G8V star t Ceti, which has a spectral 
type close to the G9V of CoRoT-7, the convective blueshifts 
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should be similar to those of the Sun, so we adopt solar val- 
ues in our simulations. In active regions, vertical convective mo- 
tions are quenched, so we observe an apparent redshift of the 
spectral lines in spotted a nd facular areas in co mparison to the 
unperturbed photosphere. iMeunier et al.l d2010t) have quantified 
this effect in the case of the Sun and we adopt their results, con- 
sidering an apparent redshift AVf = 200 m s in faculae and 
AV S = 300 m s -1 in cool spots. 

In principle, the integrated effect of convective redshifts can 
be measured in a star by comparing RV measurements obtained 
with two different line masks, one including the sh allow lines 
and the other the deep lines (cf. Meunier et al.l2010l) . In the case 
of CoRoT-7 , lacking such measurements, we apply the results of 
iGravl d2009l) and adopt solar-like values as the best approxima- 
tion. 

Considering solar convection as a template, intense down- 
drafts are localized in the dark lanes at the boundaries of the 
upwelling granules, but they contribute a significantly smaller 
flux because of their lower brightness and smaller area. While 
a consideration of those downdrafts is needed to simulate the 
shapes of line bisectors, it is beyond our simplified approach 
that assumes that the whole profile of our template line forms 
at the same depth inside the photosphere. Therefore, we restrict 
our model to the mean apparent RV variations, neglecting the as- 
sociated variations of the bisector shape and do not include the 
effect of convective downdrafts as well as those of other surface 
flows typical of solar active regions, such as the Evershed effect 
in sunspots (cf., e.g., Me unier et alJl2010h . 

The local central wavelength Aq of the k-th surface element 
at time tj is given by: A^j - /1r(1 + Vjy/c), where Vkj is its radial 
velocity: 



v kj 



-(v sin sin k sin[4 + Q,(tj - to)] + AV CS , 



(7) 



with c the speed of light, v sin i the stellar projected rotational 
velocity, Ar the rest wavelength of the spectral line, and AV CS the 
apparent convective redshift arising from the reduction of con- 
vective blueshifts in spots and faculae, which is parameterized 



as: 



AV CS = (AV S + QAVdJwh 



(8) 



where fa is the spot filling factor of the element, Q the facular- 
to-spotted area ratio, and /^y is given by Eq. (fJJ. 

We integrate the line specific intensity at a given wavelength 
over the disc of the star using a subdivision into 1° X 1° surface el- 
ements to obtain the flux along the line profile F(A, tj) at a given 
time tj. To find the apparent stellar RV, we can fit a Gaussian to 
the line profile F(A, tj), or we can determine the centroid of the 
profile as: 



A ce (tj) = 



jA'R(A,t j )dA 
J ( R(A,t j )dA ' 



(9) 



where 1 ?? = [1 -F(A,tj)/F c (tj)] is the residual profile of the spec- 
tral line computed from the ratio of the flux in the line F{A, tj) to 
the flux F c (tj) in the adjacent continuum at any given time tj. 

A single line profile computed with the above model can 
be regarded as a cross-correlation function (hereafter CCF) ob- 
tained by cross-correlating the whole stellar spectrum with a line 
mask consisting of Dirac delta functions giving t he rest wave- 
length and depth of each individual line (e.g., lOueloz et al.l 
2001). Therefore, to derive the RV from a single synthetic line 
profile is equivalent to measure the RV from a CCF, either by 
fitting it with a Gaussian, or by computing its centroid wave- 
length. A better approach would be that of simulating the whole 



stellar optical spectrum or an extended section of it in order to 
account for the wavelength dependence of the spot and facular 
contrasts as well a s of the convecti ve inhomoge neities (cf., e.g ., 
iDesort et alJl2007t iLagrange et alJl2010t IMeunier efall 12010) . 
Again, in view of our limited knowledge of the distribution of 
active regions on the stellar surface, our simplified approach is 
adequate to estimate the RV perturbations induced by magnetic 
fields. 



3.2.1. Some illustrative examples 

To illustrate how active regions affect the measurement of stel- 
lar RV, we consider the simple case of a single active region 
on a slowly rotating st ar, assuming stellar para met ers similar to 
the ca ses presented by Saar & Donahue] d 19971) and lDesort et al.l 
(120071) for comparison purposes. Specifically, we adopt v sin i = 
6.5 km s and a macroturbulence (q = 5.5 km s -1 to make 
macroturbulence effects clearly visible. For the moment, we 
do not include convective redshifts in spots and faculae because 
they were not taken into account in those simulations. We con- 
sider a single active region of square shape 18° x 18° located 
at the equator of a star observed equator-on, to maximize the 
RV variation. The spot filling factor is assumed to be / = 0.99 
and their contrast c s = 0.01, i.e., the spot effective temperature 
is much lower than the photospheric temperature. When facu- 
lae are included, we adopt Cf = 0.115 and Q — 9 to make their 
contribution clearly evident. In Fig. Q] we show the apparent RV 
and the relative flux variations produced by the active region vs. 
the rotation phase. The RV variations are computed by f itting a 
Gaussian to the line profile, as in lSaar & Donahuad 19971) . while 
by applying the centroid method we find values which differ by 
at most 10-15 percent. The greatest apparent variation is ob- 
tained in the case of a purely dark spot without any bright facula 
or reduced macroturbulence. In this case, the line profile inte- 
grated over the stellar disc shows a bump at the RV of the spot 
because the intensi ty of the local contin uum is reduced in the 
spot (see Fig. 1 in IVogt & Penrodl[l983l for an illustration of 
the effect). A comparison with the amplitudes calculated with 
the formulae of ISaar & Donahud dl997 ) or Des ort et all ([2007) 
shows agreement within ~ 10 percent. Note that the RV pertur- 
bation is positive and the flux is decreasing when the spot is on 
the approaching half of the disc. When the spot transits across 
the central meridian, the RV perturbation becomes zero and the 
flux reaches a relative minimum. Finally, the RV perturbation 
becomes negative when the spot is on the receding half of the 
disc and the flux is increasing. 

The inclusion of the radial-tangential macroturbulence pro- 
duces a decrease of the RV amplitude because macroturbulence 
is reduced in spots so that the local line profile becomes nar- 
rower an d deeper w hile its equivalent width is not significantly 
affected dGravll 988). This reduces the height of the bump in the 
integrated profile acting in the opposite direction to the perturba- 
tion produced by a dark spot, as it is evident in the case in which 
only a macroturbulence reduction is included in the model. Note 
that the wide-band flux variations are not affected by the inclu- 
sion of the macroturbulence in any case. 

The effect of the faculae is to increase the local continuum 
intensity when an active region is close to the limb, leading to a 
greater local absorption which produces a dip on the integrated 
line profile. Since the contrast of the faculae is greatly reduced 
when they move close to the centre of the disc, this has the effect 
of narrowing the RV peak produced by an active region when 
dark spots and faculae are simultaneously present, i.e., faculae 
increase the frequency of the RV variation. Specifically, the fre- 
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quency of the RV variation in the case of a single active region 
is twice the rotation frequency when dark spots or a reduction of 
the macro turbulence are present (cf. Fig. [1] upper panel, where 
we see a complete RV oscillation along one transit of the active 
region that lasts half a stellar rotation). The inclusion of solar- 
like faculae can increase the frequency of the variation up to 
four times the rotation frequency when facular effects become 
important (i.e., CfQ > 1). This is illustrated in the upper panel of 
Fig. [2] where the dotted line displays the RV variation produced 
by an active region with a sizeable facular contribution (see be- 
low) showing two complete oscillations along one transit. 

Next we consider the apparent redshifts owing to the reduc- 
tion of convective blueshifts in spots and faculae. To make the 
effect clearly evident, now we adopt a spot contrast c s = 0.665, 
a facular contrast C{ - 0.115, and Q - 4.5. For the convective 
redshifts, we assume: AV S = 500 m s and AVf = 300 ms^'.A 
smaller spot contrast, e.g., c s = 0.01, as in the previous simula- 
tions, makes the contribution of the local redshifted line profile 
coming from the spotted area very small, reducing the corre- 
sponding RV perturbation. In Fig. [2] we plot in the upper panel 
the RV variation produced by the single active region previously 
considered, with all the other parameters kept at the previous 
values (solid line). For the purpose of comparison, we also plot 
the RV perturbation obtained without convective redshifts (dot- 
ted line). In this case, the effects of cool spots and bright faculae 
almost balance each other because of the higher spot contrast, 
and the oscillation of the RV shows a high frequency. On the 
other hand, the inclusion of convective redshifts brings the fre- 
quency of the RV variation close to the rotation frequency, with 
a single maximum per rotation. This happens because the per- 
turbation is always positive (cf. middle panel in Fig. [2]), while all 
the other effects change their sign when the active region tran- 
sits from the approaching to the receding half of the stellar disc. 
Moreover, the convective redshifts do not depend per se on the 
contrast factors c s and C{ and the shift due to the faculae is ampli- 
fied by the facular-to-spotted area ratio Q. The modulation of the 
continuum flux is of course not affected by the convective shifts, 
so we obtain the same light curve independent of the value of 
AV S and AVf (cf. Fig. [2] lower panel). 

Therefore, particularly for not too dark spots (i.e., c s <: 0.4) 
and small v sin i, convective redshifts cannot be neglected in the 
simulation of the activit y-induced RV pertur bations in solar-like 
stars, as pointed out bv iMeunier et al.l (2010). In the case of the 
Sun, their inclusion not only increases the RV perturbation up to 
an order of magnitude, but also brings its dominant period close 
to the solar rotation period, while without their inclusion they 
find that the dominant period corresponds to the first harmonic 
of the solar rotation period. 

We note that the amplitude of the RV variation produced by 
an active region depends on several parameters that are poorly 
known, i.e., the latitude of the active region, the spot and facular 
contrasts, the convective redshifts, and the macroturbulence pa- 
rameter which is difficult to disentangle from rotational broaden- 
ing in a slowly rotating star such as CoRoT-7 (Leger et al. 2009). 
Moreo ver, the spot and f acular contrasts depend on the wave- 
length dLanza e t al. 2004), leading to a difference of =» 10 per- 
cent in the RV variations as derived from diffe rent orders of an 
echelle spectrum (cf., e.g., iDesort et al.ll20071) and the convec- 
tive redshift depends on the formation depth of a spectral line. 

The simultaneous presence of several active regions gives 
rise to a complex line profile distortion in the case of a slowly 
rotating star because the perturbations of different regions over- 
lap in the line profile due to the weak rotational broadening. This 
implies an additional 15-20 percent uncertainty in the determi- 
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Fig. 1. Upper panel: Apparent RV variation (upper panel) and 
relative stellar flux (lower panel) vs. the rotation phase produced 
by a single active region on a star seen equator-on. Different 
linestyles indicate different kinds of inhomogeneities inside the 
active region. Dashed: reduced macroturbulence without bright- 
ness perturbation; solid: dark spot and reduced macroturbulence; 
dotted: dark spot, bright facula, and reduced macroturbulence. 
The flux is measured in units of the unperturbed stellar flux. Note 
that macroturbulence perturbation alone does not affect the vari- 
ation of the stellar flux. 



nation of the central wavelength of the profile by the Gaussian fit 
or the centroid method, as found by comparing the RV determi- 
nations obtained with the two methods in the case of line profiles 
simulated with infinite signal-to-noise ratio and perfectly regular 
sampling. In consideration of all these limitations, the absolute 
values of the RV variations computed with our model are un- 
certain by 20 - 30 percent in the case of complex distributions 
of active regions, such as those derived by our spot modelling 
technique as applied to CoRoT-7 light curves. 

3.2.2. Radial velocity variations from spot modelling 

We can use the distribution of active regions as derived from our 
spot modelling to synthesize the corresponding RV variations 
according to the approach outlined in Sect. 13.21 Since our spot 
model assumes that active regions are stable for the time interval 
of each fitted time series Aff, the distribution of surface inho- 
mogeneities can be used to synthesize the RV variations hav- 
ing a timescale of Aff or longer. Active regions with a shorter 
lifetime produce a photometric modulation which appears in the 
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Fig. 2. Same as Fig. Q]when including (solid line) and not in- 
cluding (dotted line) the convective redshifts in cool spots and 
bright faculae in the simulation (see the text for the active region 
parameters). Note that the redshifts do not affect the stellar flux, 
so both cases coincide in the lower panel. 



residuals of the best fit to the light curve. As we shall see in 
Sect. 15. II most of the short-term variability occurs on time scales 
of 4 - 5 days, i.e., significantly shorter than the rotation period of 
CoRoT-7, so we can neglect, as a first approximation, the vari- 
ation of the disc position of those active regions due to stellar 
rotation and estimate their area as if they were located at the cen- 
tre of the disc. Specifically, first we express the residuals as the 
relative deviation AF/Fq between the observed flux and its best 
fit measured in units of the reference unspotted flux Fq; then we 
subtract from AF/Fq its mean value yu res = (AF/Fq) which cor- 
responds to a uniformly distributed population of active regions 
that do not produce any RV variation. At each given time f, we 
adopt the relative deviation \AF/Fq - fi KS \ as a measure of the 
filling factor of the active regions producing the short-term RV 
variations. In doing so, we neglect limb-darkening effects and as- 
sume that those active regions consist of completely dark spots 
(c s = 0). Finally, we compute the RV perturbation d ue to such 
brigh tness inhomogeneities by means of Eq. (1) of iDesort et alj 
(2007) obtaining: 



AVRi = 800(— -// re8 )(vsin0 
F 



l.i 



(10) 



where the RV perturbation is measured in m s 4 and the vsin/ 
of the star in km s" 1 ; for simplicity, Eq. (TTOb assumes that the 
perturbation has the same sign as AF/Fq - fi KS which should 



produce statistically correct results in the case of a sufficiently 
long time series, i.e., with equally probable positive and neg- 
ative deviations. In our case, Eq. ( TTOb gives an upper limit for 
the amplitude of the RV perturbation because it assumes that an 
active region spans the diameter of the stellar disc. This is appro- 
priate because we want to assess the impact of short-lived active 
regions on the detection of exoplanets. On the other hand, the 
effects of convective redshifts in short-lived active regions is es- 
timated as: 



AV Rc = (AV S + QAV f )(^- - tes ), 
F 



(ID 



where the sign of AVr c has been chosen to maximize the total 
short-term RV perturbation when added to that of Eq. ( fTOb . By 
adding both kinds of short-term perturbations to the long-term 
RV modulation produced by the active regions mapped by the 
spot modelling, we compute the total synthetic RV variation in- 
duced by stellar activity. 



4. Model parameters 



The fu ndamental stellar parameters are taken from iLeger et ail 
(2009) and are listed in Table Q] The limb-darkening par ameters 
a p , b p , and c p (cf. Sect. 13.11 ) have been derived from iKuruczl 
(2000) model atmospheres for r e ff = 5275 K, logg = 4.50 (cm 
s~ 2 ) and solar abundances, b y adopting the C oRoT white-band 
transmission profile given by lAuvergne et alj d2009t) . Recently, 
iBruntt et aD ((2010) have slightly revised stellar parameters, but 
adopting their values does not produce any significant modifi- 
cation as far as our modelling is concerned. The rotation period 
adopted for our spot modelling has been derived from a peri- 
odogram analysis of the light curve giving P I0t = 23.64 + 3.62 
days. The uncertainty of this period determination is derived 
from the total extension of the time series and represents an 
upper limit. As we shall see below, our spot model shows that 
the starspots have a remarkable differential rotation which con- 
tributes to an uncertainty of the stellar rotation period of ~ 6 
percent, i.e., of ±1.42 days (cf. Sect. l5.2l i. 

The polar flattening of the star due to the centrifugal potential 
is computed in the Roche approximation with a rotation period 
of 23.64 days. The relative difference between the equatorial and 
the polar radii e is 8.54 x 10~ 6 which induces a completely neg- 
ligible relative flux variation of w 10~ 7 for a spot coverage of 
~ 2 percent as a consequence of the gravity darkening of the 
equatorial region of the star. 

The inclination of the stellar rotation axis is impossible to 
constrain through the observation of the Rossiter-McLaughlin 
effect due to the very small planet relative radius R p /R s tar = 
0.0187 ± 3.0 x 10~ 4 , the small rotational broadening of the star 
vsin/ < 3.5 km s " 1 , and its intrinsic line profile variations due 
to stellar activity (lOueloz et al.ll2009l) . Nevertheless, we assume 
that the stellar rotation axis is normal to the orbital plane of the 
transiting planet, i.e., with an inclin ation of 80 ° 1 ± °3 from the 
line of sight (cf. ILeger et al .1120091) . 

The maximum time interval that our model can accurately 
fit with a fixed distribution of active regions Aff has been deter- 
mined by dividing the total interval, T = 130.8152 days, into 
Nf equal segments, i.e., Aff = T/Nf, and looking for the mini- 
mum value of Nf that allows us a good fit of the light curve, as 
measured by the^- 2 statistics. We found that for Nf < 7 the qual- 
ity of the best fit degrades significantly with respect to Nf > 7, 
owing to a substantial evolution of the pattern of surface bright- 
ness inhomogeneities. Therefore, we adopt Aff = 18.688 days 
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Table 1. Parameters adopted for the light curve and RV mod- 
elling of CoRoT-7. 



Parameter 



Ref/ 



Star Mass (M Q ) 
Star Radius (R Q ) 
T ee (K) 
log g (cms~ 2 ) 

h 

Pmt (days) 

6 

Inclination (deg) 
c s 

Cf 

Q 
Aft (d) 

vsin/ (km s _1 ) 
A/1 D (km s" 1 ) 
<ro(kms-') 
AV s (ms-») 
AV f (m s- 1 ) 



0.93 

0.87 

5275 

4.50 

0.284 

0.961 

-0.254 

23.64 

8.54 x 10 6 

80.1 

0.577 

0.115 

1.0 

18.688, 14.532 

1.7 

2.33 

1.2 

300 

200 



L09 

L09 

L09 

L09 

LalO 

LalO 

LalO 

LalO 

LalO 

L09 

LalO 

L04 

LalO 

LalO 

BrlO 

LalO 

BrlO 

G09, MelO 

G09, MelO 



" Reference s ; Brl O: iBruntt et al.1 (12 010); G09: iGravl J2009h : L04: 
lLanza etal.1 J2004h: L09: LeseretalJ d2009h : LalO: present study; 
MelO: lMeunieretal.U2O10l) . 



as the maximum time interval to be fitted with a fixed distri- 
bution of surface active regions in order to estimate the best 
va lue of the paramete r Q (see below). This confirms the result 
of lOueloz et al.l (120091) who estimated a starspot coherence time 
of ~ 18 days. However, for the spot modelling in Sects. I5TD and 
15.21 we shall adopt Nf = 9, corresponding to Atf = 14.532 days, 
which provides a better time resolution to study the evolution of 
the spot pattern during the intervals with faster modifications. 

To compute the spot contrast, we adopt the same mean tem- 
perature difference as derived for sunspot groups f rom their 
bolometric contrast, i.e., 560 K (Chap man et al.|[l994l) . In other 
words, we assume a spot effective temperature of 4715 K, 
yielding a contra st c s = 0.577 in the CoRoT passband (cf. 
lLanza et al.l2007l) . A different spot contrast changes the absolute 
spot coverages, but does not significantly affect thei r longitudes 
and th eir time evolution, as discussed in detail bv lLanza et al.l 
(2009a|). T he facular co ntrast is assumed to be solar-like with 
c f = 0.115 (lLanza et al.ll2004l) . 

The best value of the area ratio Q of the faculae to the spots 
in each active region has been estimated by means of the model 
of lLanza et all (I2003L cf. Sect. |3J)- In Fig. [3] we plot the ratio 
X 2 lx 2 ■ °f the total x 2 of the composite best fit of the entire time 
series to its minimum value y 1 . , versus Q, and indicate the 95 

/limn' s£ '' 

percent confidence l evel as derived from the F-statistics (e.g., 
lLampton et~allll976l) . Choosing Atf - 18.688 days, we fit the 
rotational modulation of the active regions for the longest time 
interval during which they remain stable, modelling both the flux 
increase due to the facular component when an active region is 
close to the limb as well as the flux decrease due to the dark 
spots when the same region transits across the central meridian 
of the disc. In such a way, a measure of the relative facular and 
spot contributions can be obtained leading to a fairly accurate 
estimate of Q. 



The best value of Q turns out to be Q - 1.0, with an ac- 
ceptable range extending from ~ to ~ 4. There is also a small 
interval of formally acceptable values between 6.0 and 6.5, but 
we regard it as a spurious outcome of the x 2 statistical fluctu- 
ations. Although a sizeable facular contribution cannot be ex- 
cluded on the basis of the photometric best fit, we shall find in 
Sect. l5.3l that Q < 4 is required to reproduce the amplitude of the 
observ ed RV variations. In the case of the Sun, the best value of 
Q is 9 ( lLanza et al.l l2007). Thus our result indicates a lower rela- 
tive contribution of the faculae to the light variation of CoRoT-7 
than in the solar case. The amplitude of the rotational modula- 
tion of the star was ~ 0.0185 mag during CoRoT observations 
and ~ 0.03 mag during the campaign organized by Oueloz et al.l 
(2009), i.e., ~ 6 - 10 times that of the Sun at the maximum of 
the eleven-year cycle. This indicates that CoRoT-7 is more ac- 
tive than the Sun, which may account for the reduced facular 
contrib utio n to its light v ariations, as suggested by Radick et al.l 
(1998) and Lockwoodet al. (2007). The ground-based photom- 
etry of lOueloz et al.l (120091) also indicates that cool spots dom- 
inate the optical variability of the star since it becomes redder 
when it is fainter. The use of the chromatic information of the 
CoRoT light curves to estimate the spot and facular contrasts 
and filling factors is made impossible by our ignorance of the 
unperturbed stellar flux levels in the different colour channels 
which are needed to disentangle the flux perturbations due to 
spots and faculae, respectively. The continuous variations of the 
observed fluxes do not allow us to fix such reference levels so 
that we cannot unambiguously attribute a given flux modulation 
to cool spots or bright faculae. Moreover, the long-term drifts of 
the fluxes in the individual colour channels complicate the esti- 
mate of the flux variations in each of them, making this approach 
unfeasible. 

To compute the RV variations induced by surface inhomo- 
geneities, we assume a line rest wavelength of 600 nm and a lo- 
cal thermal plus microturbulence broadening A/Id = 2.3 km s . 
The vsini = 1.7 + 0.2 km s is estimated fro m the stellar 
radius, inclination and rotation period following Bruntt et al. 
(|2010|) . The radial-tangential macroturbulence velocity is as- 
sumed & = 1 -2 ± 0.8 km s _I after Bruntt et al. (2010). Note that 
it is very difficult to obtain a good estimate of such a parameter 
because macroturbulence and rotational broadening are largely 
degenerate owing to the slow rotation of the star (cf . iLeger et aTl 
2009; Bruntt et al. 2010). The convective redshifts in cool spots 
and faculae, AV S and AVf, have been assumed to be similar to 
those of the Sun, following the discussion in Sect. 13.21 



5. Results 

5. 1. Light curve model 

We applied the model of Sect. l3.1l to the out-of-transit CoRoT-7 
light curve, considering nine intervals of duration Atf - 14.532 
days. The best fit without regularization (77 = 0) has a mean 
/i res = 3.862 x 10~ 6 and a standard deviation of the residu- 
als <x = 1.444 x 10~ 4 in relative flux units. The Lagrangian 
multiplier 77 is iteratively adjusted until the mean of the resid- 
uals /j res = -1.057 x 10~ 5 = -cro/ V/V, where N = 196 is the 
mean number of points in each fitted light curve interval Aff. The 
standard deviation of the residuals of the regularized best fit is 
cr = 1.951 x 10~ 4 . The composite best fit to the entire light curve 
is shown in the upper panel of Fig. |4] while the residuals are 
plotted in the lower panel. They show oscillations with a typical 
timescale of ~ 4-5 days that can be related to the rise and decay 
of small active regions that cover w 0. 1 percent of the stellar disc, 
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Fig. 3. The ratio of the x 1 of the composite best fit of the entire 
time series to its minimum value vs. the parameter Q, i.e., the ra- 
tio of the facular area to the cool spot area in active regions. The 
horizontal dashed line indicates the 95 percent confidence level 
foTx 2 /x 2 > determining the interval of acceptable Q values. 



i.e., comparable to a large sunspot group. These small active re- 
gions cannot be modelled by our approach because they do not 
produce a significant rotational flux modulation during the 23.6 
days of the stellar rotation period as they move across the disc by 
only « 60° - 75° in longitude. By decreasing the degree of regu- 
larization, i.e., the value of 77, we can marginally improve the best 
fit, but at the cost of introducing several small active regions that 
wax and wane from one time interval to the next and are badly 
constrained by the rotational modulation. Nevertheless, the os- 
cillations of the residuals do not disappear completely even for 
77 = 0, indicating that CoRoT-7 has a population of short-lived 
active regions with typical lifetimes of 4 - 5 days. 

A periodogram of the residuals is plotted in Fig. [5] and its 
main peaks correspond to periods of 5. 12 and 3.78 days, respec- 
tively, the latter close to the orbital period of CoRoT-7c of 3.69 
days. Fitting a sinusoid with this period, we find a semiamplitude 
of the light variation of 6 .47 x 1 -5 in relative flux units . The cor- 
responding RV perturbation, estimated with Eqs. ( fTUb and (fTTT i. 
has a semiamplitude of ~ 0.4 m s , an order of magnitude 
smaller than the oscillation attributed to the planet CoRoT-7c 
whose semiamplitude is 4.0 ± 0.5 m s . Concerning CoRoT-7b, 
the residual oscillations at its orbital frequency are practically 



5.2. Longitude distribution of active regions and stellar 
differential rotation 

The distributions of the spotted area vs. longitude are plot- 
ted in Fig. [6] for the nine mean epochs of our individual sub- 
sets adopting a rotation period of 23.64 days. Longitude zero 
corresponds to the sub-observer point at the initial epoch, i.e., 
HJD 2454398.0719. The longitude increases in the same direc- 
tion as the stellar rotation. This is consistent with the reference 
frames adopted in our previous studies (Lanza et al. 2009a b). To 
allow a comparison of the mapped active regions with the dips 
in the light curve, note that a feature at longitude £ crosses the 
central meridian of the star at rotation phase 360° - I. 

The spot distribution and its evolution suggest three large- 
scale active regions which show different intrinsic longitudinal 
motions. Assuming that this is due to a shear by differential rota- 
tion, we perform linear fits by assuming constant migration rates 



that are found to be 0.31+0. 12 deg/day for that fitted with a dot- 
ted line; 1.20 ±0.1 9 deg/day for that fitted with a dot-dashed line; 
and 0.61 + 0.27 deg/day for that fitted with a long-dashed line. 
The relative amplitude of the surface differential rotation, as es- 
timated from the difference between the greatest and the lowest 
migration rates, is AQ/O = 0.058 ± 0.017; this is only a lower 
limit since the spot latitudes are unknown. Since CoRoT-7 is 
more active than the Sun, its active regions may cover a latitude 
range greater than in the Sun where sunspot groups are confined 
to ±40° from the equator (see, e.g., IStrassm eier 2009). For the 
Sun, one finds AO/Q =* 0.04-0.05 by considering active regions 
confined to within the sunspot belt, i.e., within ± (35° -40°) from 
the equator. Thus, CoRoT-7 may have a surface differential ro- 
tation comparable to the Sun if its active regions are also mainly 
localized at low latitudes. 

It is interesting to compare the lower limit for the surface 
sh ear of CoRoT-7 wi th the surface differential rotation estimated 
bv lCroll et all d2006l) for the K2 dwarf e Eridani. They found two 
spots rotating with periods of 11.35 and 11.55 days at latitudes 
of 20. °0 and 31°. 5, respectively, leading to a differential rotation 
amplitude of half of that of the Sun. This is consistent with a 
weak dependence of the absolute surface shear AQ. on the angu- 
lar velocity Q, because e Eri has a rotation pe riod roughly half of 
that of the Sun. Indeed. iBarnes et alJ d2005h found AQ cc O 15 
for late-type stars at a fixed effective temperature. This suggests 
that the true pole-equator angular velocity difference in C0R0T- 
7 may be about a factor of 2 - 3 greater than observed, if the 
active regions mapped by our technique are mainly localized at 
low latitudes. 

The stable active longitudes plotted in Fig.|6]show significant 
area changes on a timescale as short as two weeks although their 
overall lifetime may exceed t he duration of the li ght curve, i.e., 
130 days. As a matter of fact. lOueloz et all (f2009) found that the 
rotational modulation of the optical flux observed during their 
ground-based campaign from December 2008 to February 2009 
matched the extrapolation of the C0R0T light curve when a ro- 
tation period of 23.64 days was assumed, although 17 complete 
rotations had elapsed from C0R0T observations. This suggests 
that the active longitudes may exist for several years, although 
some large active regions wax and wane on time scales rang- 
ing from two weeks (i.e., the time resolution of our mapping) 
to a few months. The variation of the total spotted area of our 
surface models vs. time is plotted in Fig. [7] and shows a char- 
acteristic timescale of ~ 60 - 80 days. Note that the absolute 
values of the area per longitude bin and of the total area de- 
pend on the spot and facular contrasts adopted for the modelling. 
Specifically, darker spots lead to a smaller total area while the 
effect of the facular contrast is more subtle and influences some- 
what the longitudin al distribution of the active regions (cf., e.g., 
lLanza et alj|2009al Figs. 4 and 5). 

5.3. Activity-induced radial velocity variations 

To simulate the apparent RV changes induced by the distribu- 
tion of active regions derived from our light curve modelling, 
we consider a spectral line with a rest wavelength of 600 nm, 
i.e., close to the isophotal wavelength of C0R0T observations for 
which our contrast coefficients are given. The spectral resolution 
of the line profile is A/AA = 230, 000, and 50 profiles are com- 
puted per stellar rotation to warrant a good phase sampling. The 
RV is derived by computing the centroid of each profile which 
turns out to be more stable than fitting a Gaussian which some- 
times leads to poor results in the presence of complex line profile 
distortions. First we consider the RV perturbation produced by 
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Fig. 4. Upper panel: The out-of-transit light curve of CoRoT-7 (dots) and its ME-regularized best fit for a facular-to-spotted area 
ratio of Q - 1.0 (solid line). The flux is normalized to the maximum observed flux. The vertical dashed lines mark the individually 
fitted intervals of Atf = 14.532 days. Lower panel: The corresponding residuals. 
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Fig. 5. Lomb-Scargle periodogram of the residual time series of 
Fig. 3] The dashed horizontal line marks a false-alarm probabil- 
ity of 10~ 6 . The vertical dot-dashed lines mark the orbital fre- 
quencies of the planets CoRoT-7b (right, P or b = 0.8536 days) 
and CoRoT-7c (middle, P 01 -b = 3.698 days) as well as the stellar 
rotation frequency (left, P rot = 23.64 days). 



active regions with an evolutionary timescale equal to or greater 
than Atf = 14.532 days as mapped by means of our spot mod- 
elling approach for Q — 1 and Q — (cool spots only). We plot 
the simulated RV variations vs. time in the upper panel of Fig. [8] 
(solid line for Q = 1, dotted line for 2 = 0) that shows a peak- 
to-peak amplitude up to ~ 24.5 m s -1 for Q - 1 and to ~ 20.5 m 
S for Q — 0. Small discontinuities are seen at the endpoints of 
each Aft interval owing to the change of the spot configuration 
between successive light curve fits. Their impact on the simu- 
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Fig. 6. Distributions of the spotted area vs. longitude during the 
time intervals centred on the labelled times (HJD - 2450000.0), 
adopting Q - 1.0. The plots have been vertically shifted to show 
the apparent migration of individual active regions and repeated 
beyond 0° and 360°. The dotted, dot-dashed, and long-dashed 
lines trace the apparent migration of the three most conspicuous 
features, respectively (see text for details). 



lated RV time series is small and generally comparable to the 
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Fig. 7. The total spotted area as derived from our lightcurve 
model illustrated in Figs. [4] and [6] The uncertainty of the area 
has been estimated from the errors of the photometric observa- 
tions. 



measurement errors (see below), so it is not necessary to correct 
for them. 

The observations bv lOueloz et al.l d2009l) showed a peak-to- 
peak amplitude up to ~ 40 m s -1 for the 95 percent of the mea- 
surements collected during their campaign from December 2008 
to February 2009 while the photometric optical modulation was 
a factor of ~ 1 .6 greater. Since the RV perturbation scales ap- 
proximately linearly with the spot filling factor, our computed 
amplitude for Q — 1 agrees well with the observed one, within 
» 2 - 3 percent, which is a good level of precision for this kind 
of sim ulations (cf., e.g., ISaar & Donahue! fl997t iDesort et al.l 
2007). Note that most of the active regions mapped by our spot 
modelling technique fall within ±20° from the equator, owing 
to the preference of ME regularization to put spots close to the 
sub-observer latitude to minimize their area. 

The lower panel of Fig. [8] shows the simultaneous CoRoT 
white band flux variation. In the present case, facular brighten- 
ing and macroturbulence reduction play only a minor role in the 
wide-band flux and RV variations, respectively, because Q - 1 .0 
and £o = 1.2 km s < vsin/ (cf. Sect. 13. 2\ . The dominant effect 
is the reduction of convective blueshifts in the spotted and facu- 
lar areas. The d iamonds show the sa mpling intervals of the RV 
observations bv lOueloz et al.l (120091) . Since their HARPS mea- 
surements were not simultaneous with CoRoT observations, the 
former have been arbitrarily shifted by ten days with respect to 
the latter to illustrate the effects of an examplary ground-based 
time sampling. 

To simulate more realistic RV time series, we add the short- 
term RV fluctuations associated with the residuals of our com- 
posite photometric best fit, computed according to the method 
outlined in Sect. 13.2.21 To account for HARPS measurement er- 
rors, we further add a Gaussian random noise with a standard de- 
viation of 2.0 m s -1 . In such a way, we take into account both the 
RV variations produced by active regions with a lifetime com- 
parable to Aff or longer, as mapped by our spot modelling ap- 
proach, as well as the variations induced by active regions with 
a shorter lifetime whose photometric effects appear as the resid- 
uals of the best fits. 

We compute 10 5 RV time series for 2=1, differing by 
the realization of the Gaussian noise and the initial epoch which 
is chosen from a uniform random distribution to sample all the 
parts of the synthetic RV curve in Fig.[8]and the short-term fluc- 
tuations. To allow for different initial epochs, we extend the syn- 
thetic time series by mirror reflections at their endpoints. We 
choose to compute 10 5 simulations to have enough statistics to 



sample the distribution of the RV amplitudes down to a proba- 
bility of the order of 10~ 4 . 

For each time series, we compute the Lomb periodogram 
which gives the squared amplitude of a best fitted sinusoidal sig- 
nal as a func tion of the frequency for an uneven time sampling 
(lLombll 19761) . On the other hand, the amplitudes obtained with 
the Scargle periodogram are normalized in a different way, as re- 
quired for the computation of the false-alarm probability, mak- 
ing the compar ison of the po wer level at different frequencies 
less immediate dScargldll982l) . 

The power spectrum obtained by averaging the Lomb peri- 
odograms of our 10 5 synthesized time series is plotted in the 
upper panel of Fig.[9](solid line). Strong aliases at 1.0 day 1 are 
clearly evident in the power spectrum and are due to the daily ca- 
dence of ground-based observations. The most prominent peaks 
fall at the rotation frequency and its first harmonic, as expected 
when a few dark spots dominate the RV variations (cf. Sect. l3.2.T1 
and Figs. Q] and |2). This is mainly due to the convective redshifts 
in active regions because RV time series computed without their 
effects have a mean amplitude of only ~ 18 m s and an aver- 
age power spectrum peaking at the first harmonic with a power 
~ 2.5 times lower at the rotation frequency. This happens be- 
cause all the other perturbations change sign when an active re- 
gion transits from the approaching to the receding part of the 
stellar disc, while convective redshifts are always positive. Note 
that at the orbital frequency of CoRoT-7b and CoRoT-7c, the 
power is smaller by factors of ~ 3 and ~ 4 relative to the rota- 
tion frequency, respectively. 

Since the HARPS data are not simultaneous with CoRoT 
data, we cannot compare this power spectrum with the peri- 
odogram of the HARPS RV time series. Nevertheless, i t is of 
interest to apply the technique proposed by Oueloz et al. (2009, 
Sect. 4) to filter out RV variations due to stellar activity to see 
to what extent it is capable of removing the simulated variations 
associated with stellar surface inhomogeneities. 

5.3.1 . Filtering out activity-related variations from the 
synthetic radial velocity curve 

To fit activity-induced RV variations, we model the radial veloc- 
ity v-gitj) at a given time tj as a linear combination of sinusoids 
at the rotation frequency Q. = 2n/P rot , its first and second har- 
monics: 



vrOj) = ao + y b m cos(mQtj) + c m sin(mQf ; ), 



(12) 



where the free coefficients ao, b m , and c m are determined by 
minimizing the chi square between the model and the observa- 
tions. Specifically, we apply the normal equation method and 
use LU decomposition and backsubstitution to solve the cor- 
responding setof linear equations f or the model coeffic ients 
(lPressetal.ll 1991 Sect. 15.4). As in lOueloz et all (120091) . we 
adopt P rot - 24.0 days (although the simulated time series 
were computed with a model star rotation period of 23.64 days) 
and apply the three-harmonic fit to subintervals of 20 days. 
Following Ou eloz et al] (2009), for each point of the time series 
several best fits are obtained by scanning the 20-day window 
through the full data range. These fits are then averaged to com- 
pute the best model for that point. Finally, the best fit time series 
is subtracted from the initial time series to obtain the filtered se- 
ries. A Lomb periodogram is computed for each thus filtered se- 
ries. Fig.|9]shows the power spectrum obtained by averaging the 
periodograms of all the 10 5 residual time series (dashed line). 
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Fig. 8. Upper panel: The synthesized RV variation due to stellar active regions vs. time as derived from our ME spot modelling with 
Q — 1 (solid line) and Q — (dotted line). To illustrate the effect of an exemplary ground-based sampling, open diamonds mark 
the times of HARPS observations by Oueloz et alj d2009t) after shifting their initial epoch to ten days after the beginning of CoRoT 
observations. Lower panel: The synthesized relative flux variation in the CoRoT white passband vs. time (solid line) with HARPS 
time sampling marked by open diamonds. Of course, the synthesized fluxes are the same both for models with (2=1) and without 
(Q - 0) faculae since they must reproduce the rotational modulation as observed by CoRoT 



The power due to stellar activity is reduced by factors of ~ 4 
and ~ 2 at the orbital frequencies of CoRoT-7b and CoRoT-7c, 
respectively. 

To check whether our simulated RV curve of CoRoT-7 con- 
tains spurious oscillations at the orbital frequencies of CoRoT- 
7b or CoRoT-7c, we fit a trigonometric polynomial of the form 
B sin(w ol b0 + C cos(w or bf) to each residual time series and com- 
pute the corresponding semiamplitude A(w or b) = Vfi 2 + C 2 , 
where u> 0T b is the orbital frequency of the corresponding planet. 
In the considered 10 5 simulations, A is always lower than 0.85 
m s" 1 at the frequency of CoRoT-7b and lower than 1.7 m s 
at the frequency of CoRoT-7c. The RV oscillations attributed to 
CoRoT-7b and CoRoT-7c have semiamplitudes K - 3.5 + 0.6 
and K - 4.0 ± 0.5 m s , respectively. Therefore, the probabil- 
ity that they come from the residual effects of stellar activity is 
below 10~ 4 in both cases. 

Note that, since the amplitude of the light modulation dur- 
ing the HARPS observing season was greater by a factor of 
~ 1 .6 than during CoRoT observations, the activity-induced RV 
oscillations are expected to have been greater during those ob- 
servations than those simulated from our spot modelling maps. 
Nevertheless, even increasing the amplitudes of the residual os- 
cillations by a factor of 2, the false alarm probability remains 
below 10~ 4 . 



5.3.2. Constraining the facular area 

While our photometric modelling does not exclude a facular-to- 
spotted area ratio Q > 1, we can constrain it on the basis of the 



observed RV curve. By computing spot models with Q ranging 
from 4 to 10, we find that the peak-to-peak amplitude of the RV 
variations increases from 45.7 to 120.3 m s -1 for a light curve 
relative amplitude of 0.0185 and spots located within ±20° from 
the equator. Given an observed amplitude of ~ 40 m s for a 
photometric relative modulation of ~ 0.03, we conclude that the 
amplitude of the observed RV curve can be reproduced only if 
Q < 4 which is compatible with our adoption of Q — 1 corre- 
sponding to the minimum of the x 1 of the composite photometric 
best fit. Since the RV amplitude is dominated by convective red- 
shifts in faculae and spots, it could be possible to reduce it for 
Q > 4 by reducing AV S and AVf, or decreasing the spot contrast 
c s that leads to smaller spots and hence facular areas for a fixed 
Q. However, the required decrease of AV S and AVf should be at 
least of ~ 30 percent for Q ~ 6, i.e., significantl y grea t er tha n 
allowed by the scaling relationship proposed by Gray (2009). 
Similarly, the required decrease of c s should be too large to be 
compatible with our assumption of active regions with solar-like 
contrasts. 

Our model assumes that the active regions are located close 
to the equator, according to the ME prescription that tries to min- 
imize their area as much as possible. Nevertheless, we also com- 
pute a spot model with Q - 6.5 and the constraint that active 
regions are located at high latitudes (> 60°), to study the impact 
on the RV perturbations. The resulting RV time series is very 
similar, with differences never exceeding ±6 m s , to the case 
of low-latitude spots. This happens because the reduction of the 
convective redshifts due to projection effects is compensated by 
an increase of the continuum flux perturbation due to the faculae 
since they are closer to the limb for most of the time. In sum- 
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Fig. 9. The average Lomb periodogram of 10 5 synthetic RV time series vs. frequency (solid line) and the average periodogram 
obtained after filtering the series with the three-harmonic model of Ouel oz et al.l (2009) (dashed line). The power is normalized to 
the maximum value of the unfiltered mean periodogram. The dot-dashed lines mark the rotation frequency, its first, and second 
harmonics, while the three-dot-dashed lines mark the orbital frequencies of CoRoT-7c (left) and CoRoT-7b (right), respectively. 
Note the strong aliasing associated with the daily cadence of the ground-based observations. 



mary, we conclude that a model with Q = 1 is appropriate for 
CoRoT-7, even if its active regions were not constrained to low 
latitudes. 



6. Discussion 

We analyse the white band light curve of the active late-type star 
CoRoT-7 that hosts a transiting telluric planet; a detailed analysis 
of its RV variations additionally p rovided evidence of a second 
non-t ransiting telluric planet dLeger et al.l 120091: lOueloz et al.l 
2009). Since the major source of RV variations in CoRoT-7 is 
stellar magnetic activity, it needs to be modelled in order to re- 
fine the planetary orbital parame ters and con firm the presence 
of the second planet CoRoT-7c (Ouelo zet al.ll20 09). In princi- 
ple, the RV oscillation attributed to CoRoT-7c could be spurious 
and due to the effects of stellar activity because its frequency 
is close to the fifth harmonic of the stellar rotational frequency. 
Unfortunately, the CoRoT observations are not simultaneous to 
the RV time series by Oueloz et al. (2009) so our conclusions 
based on the modelling of CoRoT data do not constrain, rigor- 
ously speaking, magnetic activity properties during the RV ob- 
servations. Nevertheless, we can derive illustrative conclusions 
from our results which can be applied to the analysis of the RV 
time series. 

Our favoured spot model has a facular-to-spotted area ratio 
Q - 1, as expected for a star with an activity level significantly 
higher than the Sun. In our simulated RV time series, most of 
the power falls at the rotation frequency and its first harmonic. 
In their analysis of the observed RV time series of CoRoT-7, 
lOueloz et al.l (12009) find that most of the power falls at the ro- 
tation frequency and at frequencies corresponding to periods of 



9.03 and 10.6 days, the latter corresponding to the first harmonic 
of the rotation frequency; another prominent peak corresponds 
to a period of 27.6 days (cf. their Table 2). The period of 27.6 
days could be associated with the rotation of spots at high lati- 
tudes given that the relative amplitude of the surface differential 
rotation of CoRoT-7 is at least of ~ 6 percent. As discussed in 
Sect. 15.21 a pole-equator angular velocity difference compara- 
ble to the Sun, i.e., « 20 percent, can indeed be present in our 
star. This is compatible with a rotation period of 23.6 days at the 
equator and of 27.6 days at high latitudes. The observed peak 
corresponding to a period of 9.03 days may therefore be the sec- 
ond harmonic of the high-latitude rotation frequency. Recently, 
Hatzes et al. (private comm.) suggested that such a modulation 
may be associated with a third planet. Our approach cannot be 
applied to support or disprove this hypothesis because our RV 
model does not include the effect of surface differential rotation 
and high-latitude active regions that may be producing the 27.6 
day periodicity and its harmonics. Therefore, we have to post- 
pone an investigation of the role of activity-induced effects for 
the 9.03 day periodicity to a future work. 

The presence of a large facular component (Q > 4) is dis- 
favoured in the case of CoRoT-7 because the amplitude of the 
RV variation, in this case dominated by convective redshifts in 
faculae, would then be too large to be compatible with the RV 
observations. 

Finally, we note that in the periodogram of the residuals 
of our ME composite best fit there is a significant peak corre- 
sponding to a period close to the orbital period of CoRoT-7c (cf. 
Fig. 0. This peak is also found in the periodogram of the time 
series of the residua ls of the best fit obtained with the method of 
lOueloz et al. (2009), so it is not an artifact associated with the 
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ME regularization. The corresponding RV semiamplitude, esti- 
mated with the method described in Sect. 13.2.21 is ~ 0.4 m s , 
i.e., about one order of magnitude smaller than the RV oscilla- 
tion produced by the planet, so these fluctuations cannot account 
for the planetary signature. On the other hand, it cannot be ex- 
cluded that those photometric oscillations, with a semiamplitude 
of only ~ 6 x 10~ 5 mag, could be induced by the interaction of 
the planet with the stellar coronal magnetic field or by a pertur- 
batio n of the stellar dyn amo by the planet, as conjectured by, 
e.g. JLanzal (120081 120091) . 

7. Conclusions 

We have applied the spot modelling method introduced by 
lLanza et al.l d2009alfbh to the lightcurve of the planet-hosting star 
CoRoT-7. Adopting a facular-to-spotted area ratio Q — 1, we 
have found evidence of three persistent active longitudes within 
which individual active regions form and decay with lifetimes 
ranging from a couple of weeks to a few months. The active lon- 
gitudes migrate in the adopted reference frame at different rates 
which suggests a surface differential rotation with a relative am- 
plitude of at least 0.058 ± 0.017. The overall spotted area reaches 
a maximum of 1 .6 percent of the whole surface if the contrasts 
of dark spots and faculae are adopted as solar-like. The residuals 
of the best fit indicate a population of small spots with areas of 
the order of 5 x 10 4 of the whole stellar surface that evolve on a 
typical timescale of ~ 4 - 5 days. Moreover, we find evidence of 
an oscillation of the optical flux with a semiamplitude ~ 6 x 10~ 5 
mag with the orbital period of the planet CoRoT-7c. 

We introduce a model to simulate the apparent RV variations 
induced by magnetic activity which includes the effect of sur- 
face brightness inhomogeneities as well as a reduction of the 
macro turbulence and convective blueshifts in active regions. We 
synthesize a time series for the RV variations of CoRoT-7 adopt- 
ing the spot model obtained from CoRoT white band light curve 
and assuming that active regions have solar-like contrasts and are 
located within ±20° from the equator. Unfortunately, HARPS 
and CoRoT observations are not simultaneous, so any model 
based on the latter does not constrain, rigorously speaking, the 
observed RV variations. Nevertheless, an extrapolation of our 
model matches the observed RV amplitude. 

We fit the RV time series with a linear combination of three 
sinusoids at the rotation frequ ency and its first two harmonics, 
following Ouel oz et alJ (I2009L Sect. 4), and confirm that their 
method is appropriate to filter out most of the activity-induced 
RV variations. The residual amplitudes at the orbital periods of 
the planets CoRoT-7b and CoRoT-7c are used to estimate false- 
alarm probabilities that the RV oscillations attributed to the plan- 
ets are spurious effects induced by stellar activity. For both plan- 
ets we find a false-alarm probability lower than 10~ 4 . 

The possible presence of a third planet with a period of 9.03 
days cannot be confirmed by means of our approach because 
this period could be associated to an harmonic of the signal with 
a period of 27.6 days which can be produced by spots at high 
latitude on a differentially rotating star. We shall consider the 
effects of surface differential rotation in future models of the RV 
variations to investigate this issue. 
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